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ABSTRACT 



We analyze the precursor of a Forbush Decrease (FD) observed with the 
Global Muon Detector Network on December 14, 2006. An intense geomagnetic 
storm is also recorded during this FD with the peak Kp index of 8+. By using 
the "two-dimensional map" of the cosmic ray intensity produced after removing 
the contribution from the diurnal anisotropy, we succeed in extracting clear sig- 
natures of the precursor. A striking feature of this event is that a weak loss-cone 
signature is first recorded more than a day prior to the Storm Sudden Commence- 
ment (SSC) onset. This suggests that the loss-cone precursor appeared only 7 
hours after the Coronal Mass Ejection (CME) eruption from the Sun, when the 
Interplanetary (IP) shock driven by the Interplanetary Coronal Mass Ejection 
(ICME) located at 0.4 AU from the Sun. We find the precursor being succes- 
sively observed with multiple detectors in the network according to the Earth's 
spin and confirmed that the precursor continuously exists in space. The long 
lead time (15.6 hours) of this precursor which is almost twice the typical value 
indicates that the IMF was more quiet in this event than a typical power spec- 
trum assumed for the IMF turbulence. The amplitude (-6.45 %) of the loss-cone 
anisotropy at the SSC onset is more than twice the FD size, indicating that the 
maximum intensity depression behind the IP shock is much larger than the FD 
size recorded at the Earth in this event. We also find the excess intensity from 
the sunward IMF direction clearly observed during ~10 hours preceding the SSC 
onset. It is shown that this excess intensity is consistent with the measurement 
of the particles accelerated by the head-on collisions with the approaching shock. 
This is the first detailed observation of the precursor due to the shock reflected 
particles with muon detectors. 
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Introduction 



A solar disturbance propagating away from the Sun affects the pre-existing population 
of galactic cosmic rays (GCRs) in a number of ways. Most well-known is the "Forbush 
Decrease" (FD), a region of suppressed cosmic-ray density located downstream of an 
Interplanetary Coronal Mass Ejection (ICME) shock. Some particles from this region of 
suppressed density leak into the upstream region and, traveling nearly at the speed of light, 



they race ahead of the approaching shock an d are observab le as a precursory 



(LC) anisotropy far into the upstream region (IBarnden 



197ll ) (INagashima et al. 



OSS cone 



19921 ). LCs 



are characterized by intensity deficits confined to a narrow pitch-angle region around the 
sunward direction along the Interplanetary Magnetic Field (IMF) and are typically visible 



4-8 h ours ahead of shock arrival associated with major geomagnetic storms (IMunakata et al. 



2000|). 



On the basis of numerical 
shock 



Leerungnavarat et al. 



simulations of the high-energy particle transport across the 



(|2003[ ) [hereafter referred to as paper 1] derived the theoretical 
constraints for the LC anisotropy. In particular, they presented the quantitative relationship 
between the angular width of the pitch-angle distribution and the interplanetary parameters 
which include the angle between the shock normal and the upstream IMF and the mean free 
path of the pitch-angle scattering of GCRs in the turbulent magnetic field. By analyzing 
the "two-dimensional map" of the cosmic ray intensity ob served with a lar g e sing le muon 



Munakata et al. 



(120051 ) reported 



hodoscope during a LC precursor period in October, 2003, 
a lead time of 4.9 hours for the LC precursor indicating a mean free path of the pitch-angle 
scattering shorter than the theoretical expectation based on the numerical simulation. 
This implies that the IMF was more turbulent in the event analyzed by them, than a 
typical power spectrum of the IMF turbulence assumed in paper 1. They also reported, 
on the other hand, a rather broad pitch-angle distribution for the event, which implies a 
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"quasi-parallel" shock according to the numerical simulation. Although this seemed to be 
conflicting with in-situ IMF and plasma data suggesting a "quasi-perpendicular" shock, it 
was discussed that this apparent conflict can be resolved by taking account of a pair of 
shocks formed by ICMEs successively ejected from the Sun. The "Alaska model" simulation 
of the event indicated that a large shock responsible for the SSC at the Earth was overtaking 
another smaller shock westward of the Earth. During the loss cone period, therefore, it 
was possible that the Earth may have been connected with the westward shock which was 
more consistent with a quasi-parallel geometry. Thus, the situation of this event was rather 
complicated introducing an additional ambiguity in the comparison of the observation with 
theoretical predictions. In this paper, we analyze another simpler LC event associated with 
a single shock. 

Cosmic ray precursors observed with a long lead time is of particular importance 
for the possible space weather forecast using cosmic ray measurements. An accurate 
measurement of such a precursor, however, requires detectors in a global network for 



continuous monitoring 



( iMunakata et al. 



2000 



over a long 


;ime 


Belov et al. 


2001) 



;ime period more than a day prior to the SSC onset 



by the Global Muon Detector Network (GMDN) (lOkazaki et all 120081 1. For accurate 
analyses of LC events, it is also necessary to properly remove the contribution from the 
diurnal anisotropy, which always exists in space with an amplitude comparable to the 



intensity deficit due to the LC anisotropy. The GMDN can also m easure precise 



diurnal anisotropy utilizing the global sky coverage of the network (lOkazaki et al. 



the 



20081 ). 



In this paper, we develop new analysis methods for removin g the influence o: 



Munakata et al. 



;he diurnal 



fcoOSh to the 



anisotropy and apply the "two-dimensional map" analysis by 
cosmic-ray precursors observed with the GMDN in December 2006 for better visualization 
of the signatures of cosmic ray precursors. 
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The outline of this paper is as follows. In §2, we briefly summarize the features of 
the event which we analyze in this paper. In §3, we apply a new analysis method to the 
data observed by the GMDN and present the best-fit analyses of a model based on the 
theoretical prediction by paper 1. Our conclusions and discussions based on the comparison 
with the theoretical prediction are given in §4. 



2. Event overview 

An X3.4 flare was first observed at 02:14 UT on December 13, 2006 and was followed 
by a halo Coronal Mass Ejection (CME) eruption shortly after at 02:54 UT. This CME was 
accompanied by a strong interplanetary (IP) shock, which traveled interplanetary space 
with an average velocity of 1160 km/s and arrived at the Ea rth causing the Storm Sudden 



Commencement (SSC) onset at 14:14 UT on December 14 (iLiu et al. 



20081 ). The shock 



speed at 1 AU was 1030 km/s as calculated from the conservation of t he plasma mass acr oss 



the s hock. It was also derived from the Rankine-Hugoniot relations ( IVinas fc Scudder 



19861 ) that the shock normal makes an angle of about 56° with the upstream magnetic 
field indicating the quasi-perpendicular shock. An intense geomagnetic storm followed the 
SSC with peak Kp index of 8+. There was no other SSC and interplanetary disturbances 
recorded at the Earth between 02:14 UT on December 13 and 14:14 UT on December 14. 
From the data recorded by the GMDN, we routinely derive the anisotropy vector (or the 
first order anisotropy) by fitting the function Iij{t) given by 

+ iy''''it){s\,^j cosuU + clj sinu;tO 

+ er^(t)c?., (1) 
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to the pressure-corrected hourly count rate If'j'^it) at universal time t in the j-th directional 
channel of the i-th muon detector in the GMDN. In this equation, Iij{t) denotes the 
contribution from the cosmic ray density /o(t) and iy^^if), i^^^if) are the best-fit 

parameters denoting three components of the anisotropy vector in the local geographical 
coordinate system (GEO), c^j j, sjjj- and c?j j are the coupling coefficients calculated by 
assuming a rigidity independent anisotropy, ti is the local time at the location of the z-th 



Okazaki et al. 



20081 ) ■ Figs.la-lb 



detector and uj = 7r/12 (for more detail of our analyses, see 
display the hourly data of the solar wind velocity and the IMF magnitude measured by 
the ACE satellite over three day period between December 13 and 15, each as a function 
of time in the day of year (DOY) on the horizontal axis. The left and right vertical solid 
lines indicate the flare onset and the SSC onset, respectively. A clear signature of the 
IP shock is seen in abrupt increases of the solar wind velocity and the IMF magnitude 
following the SSC onset. The IMF data observed during the period between the flare and 
SSC onsets show no signature of the IMF sector boundary and indicate that the Earth 
was in the "toward" sector throughout the period. Figure Ic shows hourly values of the 
derived cosmic ray density (/o(t)), whereas Figs.ld-lf display three anisotropy components 
{ix^^{'^)i ^y^^it), ^f^^(t)) transformed to the geocentric solar ecliptic (GSE) coordinate 
system. The anisotropy components shown in Figs.ld-lf are corrected for both the solar 
wind convection anisotropy and the Compton-Getting anisotropy arising from the Earth's 
30 km/s orbital motion around the Sun. In Figure Ic, a FD is clearly observed with a 
~3 % maximum depression of the cosmic ray density following the onset of the SSC as 
indicated by the vertical solid line. The average feature of the anisotropy components 
over the period preceding the SSC in Figs.ld-lf indicates ^x^^{t) < 0, C,y^^(t) > and 
C,f^^{t) > 0, which is in a qualitative agreement with the drift model prediction for the 



aniso tropy in the "toward" IMF sector during t 



fleld flOkazaki et al. 



2008 : 



Fushishita et al 



le y4 < epoch of the solar polar magnetic 



20101 ) 
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Also plotted by black points in Figure Ig is the pressure-corrected count rate recorded 
in the vertical channel of the Sao Martinho da Serra detector in Brazil (hereafter Sao 
Martinho). It is evident in this panel that Sao Martinho observed a sharp intensity decrease 
around DOY 348.4 a few hours prior to the SSC onset indicating the LC precursor. We 
note in this figure that Io{t), ^^^^it), C,y'^^(t) and ^f^^it) are all changing in response to 
the LC precursor recorded by the Sao Martinho detector. As shown later, the residual of 
the best-fitting is larger in two periods of 347.563- 347.771 and DOYs 348.271-348.563, with 
the second period containing the LC signature in Sao Martinho data in Figure Ig. This 
shows that the simple model in eq. cannot properly describe the LC anisotropy. In the 
next section, we improve this model by taking account of the contribution from the LC 
anisotropy. 

By using the ACE IMF data, we first examine the asymptotic viewing direction along 
the sunward IMF in 24 hours prior to the SSC onset in December 14, together with viewing 
directions (after correction for geomagnetic bending) available in the GMDN. Although 
there is a large fluctuation in the IMF orientation observed by ACE, we find that the 
majority of the IMF orientation is in the southern hemisphere, as expected for the sunward 
direction along the nominal Parker field in the ecliptic plane viewed from the Earth in 
December. The LC precursor, therefore, is expected to be observed in the directional 
channels viewing around the equator or a little south of the equator. This is consistent with 
the LC signature observed by Sao Martinho in the bottom panel of Figure 1. Our detailed 
analysis in the next section will show the LC precursor also observed by a small detector 
at Hobart viewing the mid-latitude in the southern hemisphere. In the following analysis, 
we use for calculating the pitch-angle of each viewing directio n in the GMD N the nominal 



19581 ). instead of 



Parker field calculated from the solar wind velocity by ACE ( iParker 
using the ACE IMF, to avoid effects of the large fluctuation in the ACE IMF orientation. 
The anisotropy of ~30 GeV GCR intensity observed with muon detectors is rather stable. 
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changing only gradually even when the ACE IMF shows a large fluctuation. 

3. Data analysis and results 

Muon detectors measure high-energy GCRs by detecting secondary muons produced 
from the hadronic interactions of primary GCRs (mostly protons) with atmospheric nuclei. 
Since relativistic muons have relatively long lifetimes (the proper half- life being 2.2 /is) and 
can reach the ground preserving the incident direction of the initiating primary particles, 
we can measure the GCR intensity in various directions with a multi-directional detector 
at a single location. The typical energy of primary GCRs modulated by FDs is ~30 GeV 
for muon detectors. The Global Muon Detector Network (GMDN) originally started with 
three detectors at Nagoya (Japan), Hobart (Austraha) and Sao Martinho da Serra (Brazil), 
each of which is multi-directional, thus allowing us to simultaneously record the muon 
intensities in various directions of viewing. These detectors have identical design, except 
for their detection area which is 36 m^ for Nagoya, 9 m^ for Hobart and 28 m^ for Sao 
Martinho, consisting of two horizontal layers of plastic scintillators, vertically separated by 
1.73 m, with an intermediate 5 cm layer of lead to absorb the soft component radiation in 
the air. Each layer consists of an array of 1 m^ unit detectors, each with a 1 mxl m plastic 
scintillator viewed by a photomultiplier tube of 12.7 cm diameter. By counting pulses of 
the 2-fold coincidences between a pair of detectors on the upper and lower layers, we can 
record the rate of muons from the corresponding incident direction. The multi-directional 
muon telescope comprises various combinations between the upper and lower detectors. 
The Field of View (FOV) of the vertical channel has a broad geometrical aperture of 
~ ±30° (±tan~^( 1.0/ 1.73)). A precise measurement of LCs using muon detectors became 
possible only recently, when the GMDN capable of continuously monitoring the sunward 
IMF direction was completed in March 2006 by adding another 9 m^ detector at Kuwait 
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University. Unlike the other three detectors, the Kuwait Muon Telescope (hereafter Kuwait) 
consists of four horizontal layers of 30 proportional counter tubes (PCTs). Each PCT is a 
5 m long cylinder with a 10 cm diameter having a 50-micron thick tungsten anode along 
the cylinder axis. A 5 cm layer of lead is installed above the detector to absorb the soft 
component radiation. The PCT axes are aligned geographic east-west (X) in the top and 
third layers and north-south (Y) in the second and bottom layers. The top and second layers 
form an upper pair, while the third and bottom layers form a lower pair. The two pairs are 
separated vertically by 80 cm. Muon recording is triggered by the fourfold coincidence of 
pulses from all layers and the incident direction is identified from X-Y locations of the upper 
and lower PCT pairs. This is approximately equivalent to recording muons with two 30x30 
square arrays of 10 cmxlO cm detectors separated vertically by 80 cm. The geometrical 
aperture of the FOV in this detector is ~ ±7° (± tan^^(0.1/0.8)) for the vertical channel. 
The muon count is recorded in each of 23x23=529 directional channels which cover 360° of 
azimuth angle and 0° to 60° zenith (for more detail of GMDN, see 



Okazaki et al 



20081) 



The number of viewing directions available from the GMDN was drastically 
increased from the conventional one by installing a new recording system using the Field 
Programmable Gate Arrays (FPGAs), with which we can count muons for all possible 



coinc idences between a pair of unit detectors on the upper and lower layers ( lYasue et al. 



20031 ). If we have a square n x n array of unit detectors aligned to the north-south 
(or east-west) direction in the i-th muon detector, it is possible to analyze the pressure 
corrected muon rates I°^^i-^{t){k^l = —n + 1, ...,0, ...,n — 1) in total (2n — 1) x (2n — 1) 
directional channels at the time t, where positive (negative) k and / represent eastern 
(western) and northern (southern) incidence, respectively, with k = I = corresponding to 
the vertical incident. In order to visualize the directional distribution of GCR intensity, we 
plot I°^^i^{t) as a function of k and / in a color-coded format, which we call the "2D map" 
(two dimensional map of the cosmic ray intensity). The number of directional channels 
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used in this paper is 25 (5x5) from Hobart, 121 (11x11) from Kuwait and 49 (7x7) from 
Sao Martinho, respectively. We cannot apply this technique to Nagoya data in this paper, 
as the FPGA recording system was installed in the Nagoya detector only in May, 2007 and 
not available for the FD event in December 2006. In this paper, therefore, we analyze only 
the conventional 1 7 direction data from Nagoya (for the conventional directional channels 



in the GMDN, see lOkazaki et al.l 



20081 ). Thus, the total number of the directional channels 



analyzed in this paper is 212 from four detector systems in the GMDN. 



3.1. Analysis method and the sky-maps of the observed anisotropy 

As stated in §1, we need to accurately remove from the data the contribution from the 
diurnal anisotropy (DA) for precise analysis of the LC precursor. This was not possible for 
us before completing the GMDN with which we can precisely measure the DA utilizing the 
global sky coverage of the network. By using the 24-hour trailing moving averages (TMAs) 
of Io{t), i^^^{t), iy^^{t) and ^f^*^(t) in eq. ([1]), we calculate the contribution from the 
DA to /i(fc,z) if) at a universal time t as 

/i(fc,oW = ^\k,i)if)+lx^° ii){(Ai(^k,i) cos uti - sin uti) 

+^y^^{t){s\i(k,i) cos uti + c\i^k,i) sinwti) 
+er^(t)c?,(,,), (2) 

where -^°(jt;)(t) is the 24-hour TMA of the contribution from the cosmic ray density (/o(^)) 
and Io{t), i^^°{t), ^y^'^it) and i^^°{t) are the 24-hour TMAs of the best-fit parameters 
in eq. ([1]) calculated as 

Jo(t)=Ej_23/o(t)/24, (3) 
er^(t)=S^_23er°(t)/24, (4) 
er^(t)=S^_,3er°(t)/24, (5) 
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er^(t)=sL23er°(^)/24. (6) 

By subtracting Ii{k,i)(t) in eq. ([2]) from the observed I^l'i^i^it), we get the directional intensity 
distribution A/°^^*^j(t) free from the DA as 

We cannot use /o(t), ^y^°{t), ^f^°{t) in eq. ([2]) instead of /o(t), 

iy^^if), i^^^it)^ respectively, because the LC signature recorded in a large detector like 
Sao Martinho makes the best-fit Ii(k,i){t) too close to the observed value and consequently 
leads to too small AJ°^;^,)(t) in eq. 0. In Figure Ig, Ii{k,i){t) and AJ°^;,)(t) for 

the vertical channel of Sao Martinho are displayed by solid circles, a gray curve and open 
circles, respectively. It is seen that the gradual variation of I°^^i-^{t) due to the DA {Ji[k,i)if)) 
is successfully removed in AI°l'^i^{t). We have calculated using 8-hour and 12- hour TMAs 
for Ii(k,i) (t) instead of 24-hour TMA and have confirmed that the results remain essentially 
unchanged. Note that A/°^^'';^(t) in eq. (JTj) is derived using the "trailing" average and is not 
affected by the variation occurring after t. This is an important issue for possible real time 
forecasting. 

To visualize the LC signature more clearly by suppressing statistical fiuctuation which 
is larger in the inclined channels, we also use, instead of A/°^^^^j(t), the "significance" 
sli{k,i)(t) defined as 

sI^^^^i){t) = AI°^l^^{t)/auk,i), (8) 

where (Ji(^k,i) is the count rate error for the {k, I) directional channel in the i-th detector. 
Figure 2a shows the 2D maps of in eq. ([8]) for Sao Martinho on December 13 

and 14 during 36 hours preceding the SSC onset at 14:14 UT on December 14 (DOY 
348.593), whereas Figure 2c shows sl°^^i-^{t) for Hobart in the same period. Each small 
square panel in these figures displays sl°^^i^^{t) observed in one hour in a color-coded format 
as a function of k and / on the horizontal and vertical axes denoting the east-west and 
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north-south inchnations of the viewing direction, respectively. Red color denotes the excess 
intensity in each {k, I) pixel relative to the ominidirectional intensity in an entire FOV, 
while blue color denotes the deficit intensity. In these figures, we set color scales ranging ± 
5 for Sao Martinho and ± 3 for Hobart taking account of larger statistical error in Hobart 
due to its smaller detection area. Also shown by the white curve in each panel is the 
contour line of the pitch-angle measured from the sunward IMF direction and calculated 
for cosmic rays incident to each [k, I) pixel with the median primary energy appropriate 
to that pixel (we assume 1/P rigidity spectrum for the LC anisotropy throughout this 
paper). It is seen that the zero pitch-angle region is first captured in the FOV by Sao 
Martinho in DOYs 347.354-347.563 (08:00-13:00 UT, December 13), and then by Hobart 
in DOYs 347.771-347.979 (18:00-23:00 UT, December 13) and by Sao Martinho again 
in DOYs 348.354-348.563 (08:00-13:00 UT, December 14) according to the Earth's spin. 
Hereafter, these three 6-hour periods are referred to as the LC periods. In Figure 2a for Sao 
Martinho, there is a clear LC precursor seen as a deficit intensity at around 0° pitch-angle 
with the minimum intensity about -0.8 % in DOYs 348.396~348.521 (09:00~12:00 UT on 
December 14) during one of the LC periods. A striking feature of this event is that a weaker 
LC signature is also seen one day earlier in DOYs 347.437~347.521 (10:00~12:00 UT on 
December 13) during another LC period. This suggests that the LC precursor already 
existed only 7 hours after the CME eruption at 02:54 UT on December 13 , when the IP 



20081). If this is the 



shock driven by an ICME located at 0.4 AU from the Sun (iLiu et al. 
case, the signature also should have been observed during a period in between these two LC 
periods by other detectors viewing the eastern sky neighboring the FOV of Sao Martinho. 
This is actually seen in sl^l'^^^it) for Hobart in Figure 2c in DOYs 347.771 ~348.021 (18:00, 
December 13~00:00 UT, December 14). In DOYs 348.188~348.271 (04:00~07:00 UT on 
December 14), a weak LC signature is also seen in Figure 2e for Kuwait viewing the sky 
in between the FOVs of Sao Martinho and Hobart. It is rather surprising that the LC 
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signature passing along the southern edge of the FOV of Kuwait can be seen even with such 
a small detector (thanks to a better directional resolution of the muon hodoscope). We add 
to note that a large muon hodoscope GRAPES-3 in operation at Ooty in southern India 
also recorded a clear LC precursor on ~09:00 LT (03:30 UT) in D ecember 14, just before 



20061 ) [Dr. H. Kojima, private 



the signature recorded in Kuwait in Figure 2e ( iNonaka et al.l 
communication]. These all give observational supports for the picture that the LC signature 
was continuously existing and was successively observed with the multiple detectors in the 
GMDN according to the Earth's spin. 

There is also a clear intensity excess (EX) seen with the maximum intensity about 
+0.5 % at 30°-90° pitch-angle particularly in DOYs 348.188-348.563 (04:00-13:00 UT on 
December 14) in Figs. 2a and 2c, partly overlapping the LC signature in Figure 2a for Sao 
Martinho. Such an intensity excess is expected from th e ground based measurement of 



GCR s reflected by the IP shock approaching the Earth (IDorman et al. 



1995 



Belov et al. 



19951 ). The numerical model for GCR transport across the shock also predicts such an 
EX anisotropy superposed on the LC distribution (paper 1). The fractional energy gain 
for a GCR traveling along the IMF with an energy E after the reflection by the IP shock 



approaching with a velocity is calculated as 



AE/E ^2Vscos9Bn/c, 



(9) 



where Obu is the angle between the shock normal and the upstream IMF and c is the speed 
of light. The intensity excess expected from this energy gain is estimated as 



AI/I = -fAE/E ^ 2-fVscos9Bn/c, 



(10) 



where 7 is the power-law index of the GCR ener gy spectrum, w hich we set to 2.7. By 
using the IP shock velocity (1030 km/s) at 1 AU (Liu et al.ii2008l ) for Vs and by tentatively 
assuming Obu = 60°(30°), we get ~+0.9 %(+1.6 %) for AI/I. 



- 15 - 



We confirm that there is no notable excess or deficit intensity in sl°^^i-^{t) except 
that in the periods mentioned above. In the following subsection, we will present detailed 
comparisons between the observation and the theoretical model for the LC and EX 
anisotropics. 



3.2. Model anisotropy and best-fitting to the observation 

Figure 3 shows the pitch-angle distribution of the GCR intensity. Open circles in this 
figure display a sample numerical pitch-angle distribution near an oblique interplanetary 
shock taken from Fig. 5 of paper 1. In this figure, the GCR intensity is plotted as a function 
of the pitch-angle cosine (//). A sharp LC around /i = 1 (zero pitch-angle) and a broad 
excess intensity around /i = 0.5 (60° pitch-angle) are both evident in this figure. We model 
this distribution by a superposition of the LC and EX anisotropics, shown by plus signs 
and crosses in Figure 3, respectively, as given by 

f{e, P, t) = CLc{t)fLc{0, P) + CExim + COS 6), (11) 

where 

Ac(«./') = (^)-'exp(-g). (12) 

and 6 is the pitch-angle measured from the sunward IMF direction, P is the rigidity of 
GCRs in GV, Cicif) (< 0) and CExif) (> 0) are amplitudes of the LC and EX anisotropics, 
respectively, and 9q is the angle parameter denoting the width of the LC anisotropy. The 
1/P dependence of the LC amphtude on P is assumed to follow the average P-dependence 
of the size of the FD. It is seen in Figure 3 that the model anisotropy shown by a solid 
curve reproduces well the numerical pitch-angle distribution by paper 1. 

By using / defined in eq. ( ITTj) . we calculate the expected intensity Il^lnit) for the (/c, /) 
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pixel in the i-th. detector, as 

+4''''''(^)(sL(fc,o cosuU + cL(fc,o sin cut,) + 
J pent Ni(^k,i){P)dP 

i{k,l) 

where Ni(^k,i){P), representing the number of muons produced by primary particles with 
rigidity P and recorded in the {k, I) pixel in the z-th detector, is cal culated by utilizing th e 



response function of muons in the atmosphere to primary particles ( iMurakami et al.lll979l ) 
and P^^^i^ represents the minimum (cut-off) rigidity of primary cosmic rays to produce 
muons recorded in the {k, I) pixel. For detectors (Nagoya, Hobart and Sao Martinho) in 
which each {k, I) pixel has a wide FOV resulting from the 2-fold coincidence between 1 
m^ unit detectors, we perform the integration in the last term for each combination of 
the virtual 0.1x0.1 m^ sub-detectors and then sum all values to get an integration. We 
perform only one integration for each {k, I) pixel in Kuwait. The angular resolution of this 
integration is, therefore, 3.3° for Nagoya, Hobart and Sao Martinho and 7.1° for Kuwait. 
Note that the EX anisotropy, which is essentially the first order anisotropy along the IMF 
as described in eq. (IIID, is included in ^^(t), ^J^^^(t) and ^^(t), while the LC 
amplitude Cicif) explicitly appears in the model intensity in eq. (fT3ll . We will derive later 
CExit)i therefore, by calculating the component of ^'^^^(t) parallel to the IMF. We repeat 
the calculation of lil^li-^it) in eq. f|T3|) changing 6^0 in every 1° step and determine the best-fit 
parameters CLcifyi Poit), ^"^^'^(t) which minimize 5" defined as 



S 



\ m=l ^ 



M N (jobs (j- \ _ Tcal (j. \\2 
m=lj,fcJ=l H«.') 



where s{tm) is the hourly residual of the best-fitting at the time and N and M are 
the total number of directional channels and hours used for the best-fit calculations, 
respectively. Note that all the best-fit parameters Cici^)-, ^oi^)^ $,'^^^{t) at each time t 
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are uniquely determined by the linear least-square method for each value of Oq. We carry 
out this best-fit calculation for a total of 36 hours (DOYs 347.104~348.563) preceding the 
SSC onset and obtain the best-fit parameters with the minimum residual of 5* = 1.22 and 
9o = 35°. As we did for deriving sl°l'^i-j(t) in eq. ([8]), we subtract Ii(k,i)(t) in eq. ([2]) from 
^i{ki)i^) S6t the reproduced intensity distribution A/|j.'^';-)(t), as 

A/;(tz)W=/l(t)W-^^(M)W- (15) 

The 2D maps {sl'r^li^it)) reproduced from the best-fit parameters are shown in Figs. 2b, 2d 
and 2f for Sao Martinho, Hobart and Kuwait, respectively. It is seen that both the LC and 
EX features are successfully reproduced with the best-fit parameters, although we don't see 
such a clear resemblance for Kuwait in Figure 2f due to the poor statistical significance of 
the data from this small muon hodoscope. 

Figure 4 displays the best-fit parameters, each as a function of time measured from 
the SSC. In Figure 4a, the hourly residual s{tm) is shown by black circles, together with 
the residual in the conventional best-fit analyses in Figure 1 shown by gray circles. It is 
seen that the best-fitting is improved by the new model in eq. f|T3|) as indicated by the 
reduced residual, particularly when the best-fitting in Figure 1 fails with large residual due 
to the LC signatures observed by the GMDN. Black circles in Figure 4b display the LC 
amplitude {CLc{t)) obtained for the LC periods when the zero pitch- angle is monitored by 
the GMDN, while gray circles display the amplitude obtained when the zero pitch-angle 
is out of the FOV of the GMDN. There is a large fluctuation seen in CLc{t) in this figure 
because we cannot make accurate best-fitting when the LC is out of the FOV. It is clear, 
however, that the LC amplitude {CLc(t)) gradually decreases toward ~ — 6 % at the SSC 
onset. This LC amplitude is almost twice the FD size of ~-3 % (see Figure Ic), indicating 
that the maximum intensity depression behind the IP shock is much larger than the FD size 
recorded at the Earth. Also shown in the remaining panels are the best-fit parameters Io(t) 
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and $,"^^'-^{t) transformed to the GSE coordinate system. In Figs.4d-4f, we plot parameters 
(A/o(t), A^^^^{t), Ae^^^(t), Aef^^(t)) respectively calculated, as 

A/o(t) = /^(t)-/o(t), (16) 
Ae^(t)=er^(t)-er^(t), (17) 
Aef^(t)=er^(t)-ef^(t), (18) 

Aef^^(t) = ef^^(t)-eT"^(t), (i9) 

after subtracting the 24 hour TMAs of the conventional best-fit parameters representing 
the contribution from the diurnal anisotropy. We note A^^^^{t) and A^^'^^(t) in Figs.4d 
and 4e showing positive and negative deviations from zero, respectively, during ~10 hours 
preceding the SSC This is consistent with the enhancement of the excess intensity from the 
sunward IMF direction, as expected from the measurement of the shock reflected particles. 
This is conflrmed in Figure 4g showing the anisotropy from the sunward IMF direction. 
This parallel anisotropy shows the maximum deviation of ~1 % at the SSC onset which 
is consistent with the expectation from the measurement of the shock reflected particles 
as discussed in eq. ( ITOl) . We flnally note the density (A/o(t)) also showing a positive 
deviation similar to the parallel anisotropy. The maximum deviation is ~1 % at hour 
from the SSC and is almost the same as the maximum deviation of the parallel anisotropy, 
giving a support to our model function of the EX anisotropy in eq. ( !TT|) . wherein the 
second term of the right hand side contains two terms independent and dependent on the 
pitch-angle (6) with the same amplitude at ^ = 0°. In the next section, we will compare 
these best-flt parameters with the theoretical predictions by paper 1 and discuss their 
physical implications. 
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4. Summary and discussions 

We have analyzed cosmic ray precursors of a Forbush Decrease (FD) observed on 
December 2006 with the GMDN monitoring the directional intensity of ~50 GeV galactic 
cosmic rays. This Forbush Decrease was caused by an IP shock associated with a CME 
which erupted from the Sun shortly after an X3.4 flare on December 13 and arrived at the 
Earth in December 14. An intense geomagnetic storm was also recorded during this FD 
with the peak Kp index of 8+. There was no other SSC and/or interplanetary disturbances 
recorded at the Earth during a period between the onsets of the flare and the SSC, which 
is analyzed in this paper. By using the 2D maps of the cosmic ray intensity produced by 
removing the contribution from the diurnal anisotropy, we found a clear signature of the 
LC anisotropy which was observed as a deficit intensity from the sunward IMF direction. 
The significant LC signature was first recorded by Hobart at ~ 20 hours before the SSC 
and then by Sao Martinho with a larger amplitude at ~ —6 hours. A weak LC signature 
was also seen in Kuwait viewing the sky in between FOVs of Hobart and Sao Martinho. 

The "half-width" opening angle 9hw of the LC is defined in paper 1 as the pitch-angle 
at which the intensity decrease (relative to the omnidirectional intensity) has reached 
half its maximum value. According to this definition, 6hw is calculated to be 36.4° using 
^0 = 35° obtained in §3.2. By linearly interpolating the numerical relationship between the 
loss-cone width {6hw) and the angle between the shock normal and the upstream magnetic 
field {Obu) given by paper 1 (see Table 2 in paper 1), we find 9hw = 36.4° corresponding to 
Obu = 23° assuming the "local" slope of the power spectrum of the IMF turbulence (g = 0.5) 
appropria t e for the muon detector. This Obu is less than half of Obu = 56° obtained by 



LiuetaL 



(120081 ) (see §2). One possible source of this discrepancy is that a large error might 
associate with the non-linear best-fit parameter 6q obtained in §3.2. If we evaluate the error 
by calculating 6q with which S"^ reaches at S*^ = 2.49 = 1.22^ + 1, that is, +1.0 more than 
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S'^ with minimum S, 9q with the error is estimated to be 9^ = SS^Jg" which corr esponds to 



Ofin = 23l'^3o. This Obti with a large error is not inconsistent with Obu = 56° by 



LiuetaL 



( 120081 ) ■ We are not sure whether this error estimation is appropriate or not, but we also 
carried out the best-fit calculation in §3.2 with 6q fixed at 23° corresponding to 6Bn = 56° 
and confirmed that the minimum residual 5* and the best-fit parameters except 6q are 
almost unchanged from those in Figure 4. The reproduced 2D maps were also very similar 
to those in Figs. 2b, 2d and 2f. This implies that Osn cannot be determined accurately by 
the present analysis method for the GMDN data. One reason for such a large error in Osn 
(or 6q) is the poor angular resolution of the observed incident direction in detectors (Hobart 
and Sao Martinho) with which the LC anisotropy was recorded. This is different from 
another LC event on October 2003 which was obser ved with a large muon hodoscope with 



the better angular resolution ( iMunakata et al 



20051 ) . It is also noted that the dependence 



of 9hw (i-G- ^o) on Obu becomes very small for Obu above ~ 30° according to paper 1. This 
means that one has to measure very accurately for precisely determining Obu particularly 
when Obu is large in case of the quasi-perpendicular shock. 

A striking feature of this event is that a weak LC signature was also recorded by Sao 
Martinho more than a day earlier on December 13 at ~ —28 hours. To examine the lead 
time of this LC precursor, we fit to Cicit) in Figure 4b an exponential function of time t 
(measured from the SSC) defined as, 

CBc{t) = CLc{0)exp{U, (20) 

Jo 

and obtain Clc{0) = —6.45 % and Tq = 15.6 hours. For this fitting, we used only Cicif) 
obtained when the zero pitch-angle is monitored by the GMDN (as displayed by solid circles 
in Figure 4b). This lead time Tq for 30 GeV particles corresponds to the "decay length" I 
of the LC, as 

/ = ToVs/cosObu ~ 0.42AU, (21) 
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with Obu = 23° and Vs = 1030 km/s flLiu et al.l 120081 ) . By using a ratio l/X = 0.16 derived 



for q = 0.5 and 6hw = 36.4° from interpolating Table 2 in paper 1, we estimate the parallel 
mean free path for the pitch-angle scattering (A) to be ~2.6 AU, which is almost two times 
longer than 1.5 AU estimated by paper 1 for the muon detector. This indicates that the 
IMF was more quiet in this event than a typical power spectrum assumed for the IMF 
turbulence by paper 1. These conclusions with q = 0.5 remain unchanged even if we choose 
q = 1.0, with which we get Osn ~ 25°, / ~ 0.43 AU and A ~ 3.2 AU. 

We have also found the excess intensity from the sunward IMF direction clearly 
observed during ~10 hours preceding the SSC We have shown that this excess intensity 
is consistent with the measurement of the particles accelerated by the head-on collisions 
with the approaching shock. This is the first detailed observation of the precursor due to 
the shock reflected particles with muon detectors. It is interesting to note that there is also 
an enhancement seen in the anisotropy perpendicular to the IMF in Figure 4h. We think 
that this perpendicular anisotropy is probably arising from the drift anisotropy expressed 
by the vector product between the IMF and the spatial gradient of the cosmic ray density, 
which directs toward the shock due to the higher population of the reflected particles at 
the shock. By using the observed perpendicular anisotropy together with the IMF data, 
therefore, we can deduce the temporal variation of the cosmic-ray density gradient and infer 
the geometry of the IP shock. In our separate papers, we actu ally succeeded in deriv i ng th e 



geometry of the magnetic flux rope in ICMEs by this method (iKuwabara et al. 



2004 . 



20091). 



The precursor observed with the long lead time, like the event analyzed in this 
paper, is of particular importance for the possible space weather forecast using cosmic 
ray measurements. For the accurate observation of such event, however, we need further 
improvement of the GMDN. First, the incomplete sky-coverage of the GMDN allowed us to 
analyze the best-fit parameters during only half a period in Figure 4, when the sunward 
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IMF direction was in the FOV of the GMDN. Second, the insufficient detection areas of 
Hobart and Kuwait increased the statistical error and introduced a non-uniformity into the 
response of the GMDN to the LC precursor. We beheve that such a non-uniform response 
also contributes to the large fluctuations in the best-flt parameters in Figure 4. We are 
now planning to overcome these technical problems by expanding the detection areas of 
Hobart and Kuwait and also by installing new detector(s) to expand the FOV of the GMDN 
preparing for the next solar maximum expected at around the year 2013. 
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Grants-in-Aid for Scientific Research from the Ministry of Education, Culture, Sports, 
Science and Technology in Japan and by the joint research programs of the Solar- Terrestrial 
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Telescope are supported by the Kuwait University grant SP03/03. We thank N. F. Ness for 
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Fig. 1. — Hourly mean solar wind and cosmic ray data observed during three days between 
13 and 15 December, 2006. Each panel from the top shows (a) the ACE-Level2 solar wind 
velocity, (b) the ^C£^Level2 IMF magnitude, (c) the best-fit GCR density (/o(t)), (d-f) 
three components {Cx^^{t), Cy^^{t), Cz^^{t)) of the GCR anisotropy observed by the GMDN 
and (g) the pressure-corrected count rate recorded in the vertical channel of the 

Sao Martinho detector in Brazil, as a function of time in the day of year (DOY) on the 
horizontal axis. The anisotropy components in this figure are transformed to the GSE 
coordinate system and corrected for both the solar wind convection anisotropy and the 
Compton- Getting anisotropy arising from the Earth's orbital motion around the Sun. The 
onset times of the X3.4 flare at 02:14 UT on December 13 and the SSC at 14:14 UT on 
December 14 are indicated by vertical sohd lines. The gray curves in c-f display the 24 
hour TMAs of data shown by solid circles (Io{t), ^^^^{t), ^^^^{t), ^f^^{t), in eqs.(3)-(6) 
). The gray curve in g shows Ii{k,i){t) for the vertical channel of Sao Martinho expected 
from the 24- hour TMAs of the best-fit parameter (see §3.1 in text), while open circles in g 
display AI°^^i-^{t) (see text). The gap in the solar wind velocity data between 05:00-18:00 UT 
on December 13 (DOYs 347.229-347.354) is filled by using the ^C£;Level2 alpha-particle 
velocity. 



-26- 



Fig. 2. — 2D significance maps observed by the GMDN on December 13 and 14 preceding 
the SSC onset. Panels a, c and e display the 2D maps isl°lj^i^{t)) observed by Sao Martinho, 
Hobart and Kuwait, respectively, while b, d and f show the maps {sl^^^i-^{t)) reproduced from 
the best-fit parameters obtained in §3.2. Each small square panel in this figure displays the 
"significance" in one hour in a color-coded format as a function of k and I on the horizontal 
and vertical indicated at the left and bottom edges of each figure. Red (blue) color 

denotes the excess (deficit) intensity in {k, I) pixel relative to the ominidirectional intensity 
in the FOV. The color scale is set ± 5 for Sao Martinho and ± 3 for Hobart and Kuwait as 
indicated by color bars at the right bottom corners of panels b and d. Also shown by the 
white curve in each square panel is the contour line of the pitch-angle measured from the 
sunward IMF direction and calculated for cosmic rays incident to each {k, I) pixel with the 
median primary energy appropriate to that pixel. All maps during 36 hours are shown in 
a-d for Sao Martinho and Hobart, while maps by Kuwait are shown only for 6 hours when 
the weak LC signature is visible with this small detector. 
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Fig. 3. — Model pitch-angle distribution fo r the best-fit analyses 



numerical GCR intensity read from FIG. 5 of 



Leerungnavarat et al. 



Open circles display the 



( 120031 ) as a function of the 



pitch-angle cosine on the horizontal axis, while the solid line shows our model distribution 
which is the sum of the Gaussian (plus signs) and pitch-angle cosine (crosses) functions 
denoting the LC and the EX anisotropics, respectively. In this figure, the model distribution 
is normalized to 1.239 at /i = —1. 
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Fig. 4. — Best-fit parameters in the model pitch-angle distribution. These parameters are 
obtained from the best- fitting to the observed 2D maps in Figure 2. Each panel from the top 
displays by solid circles (a) the hourly residual (s(tm))) (b) the amplitude of the LC anisotropy 
{CLc{t)), (c) the density (A7o(^)), (d-f) three components of the first order anisotropy in the 
GSE coordinate system (A^^^^(i), A^^'^^{t), A^^^^{t)) after subtracting 24-hour TMA, (g) 
and (h) the anisotropy component parallel and perpendicular to the sunward IMF direction, 
respectively, as a function of time measured from the SSC. For comparison, the residual 
in the conventional best-fit analysis is also plotted by gray circles in a. Solid circles in b- 
h display the parameters obtained when the sunward IMF direction is monitored by the 
GMDN, while gray circles display parameters obtained when the sunward IMF direction is 
out of the FOV of the GMDN. Each gray curve in b, c and g is the best-fit to the sohd circles 
with the exponential function of the time (see text). 
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